in 



o 



Spectra of absolute instruments 
from the WKB approximation 

Tomas Tyc 

Faculty of Science and Faculty of Informatics, Masaryk University, 

Kotlafska 2, 61137 Brno, Czech Republic 

April 3, 2013 



Abstract 



o 

pg We calculate frequency spectra of absolute optical instruments using the WKB approx- 
imation. The resulting eigenfrequencies approximate the actual values very accurately, in 

c/j some cases they even give the exact values. Our calculations confirm results obtained pre- 

^ ^ viously by a completely different method. In particular, the eigenfrequencies of absolute 

'q^ instruments form tight groups that are almost equidistantly spaced. We demonstrate our 

O method and its results on several examples. 

O 

"^ 1 Introduction 



Absolute optical instrument (AI) is a device that provides a perfectly sharp image of all points 

>■ in some spatial region [1]. The simplest AI is a plane mirror that gives a virtual image of a 

whole half-space. Another, beautiful example of an AI is Maxwell's fish eye, discovered by 

\Q J- C. Maxwell in 1854 ||2|, that images sharply the whole space, and all rays form circles. In 



recent years absolute instruments attracted an increased interest which has led to proposing 



^ new devices of various types, e.g. AIs that perform imaging of the whole space, of optically 

^r^ homogeneous regions [21 0], or provide magnified images [5]. A general method has been 

^ proposed in [4] for designing spherically symmetric AIs. This research was based on geometrical 

>■ optics. 

S^ Recently AIs attracted attention also from the point of view of wave optics. It was shown 

H both theoretically [6] and experimentally [71 IH] that these devices can provide subwavelength 

resolution, although this claim has raised controversy [9l [TOl [11] and it is still not clear to what 
extent such a super-resolution can be used practically [12]. A different question was addressed 
in [T3], namely what are the general characteristics of the spectrum of eigenfrequencies of 
AIs. It was shown that the spectrum consists of tight groups of levels with almost equidistant 
spacing between them. This finding was based on an analysis of a light pulse propagating in 
the AI and on the assumption that a short pulse emitted at some point can be absorbed at the 
image point during a short time as well. Numerically calculated spectra of various AI confirmed 
this theoretical result very well. 

In this paper, we investigate the spectra of absolute instruments by a completely different 
method. Employing the WKB approximation with the Langer modification and using one of 

1 



the general properties of AIs, we confirm in a different way the previously known results about 
their spectrum [13]. Our method has two advantages compared to the previous one: it enables 
to calculate not just the spacing of the level groups but also their offset, and it allows to treat 
the situations where a mirror is used in the device. We verify our results by comparing the 
calculated spectra with numerical values for several examples of AIs. 

The paper is organised as follows. In Sec.[2]we recall absolute instruments and discuss some 
of their properties. In Sec. [3] we employ the WKB method for calculating the spectra of radially 
symmetric media and in Sec. |4]we illustrate the results on particular examples of AIs. In Sec. [5] 
we analyse the situation in AIs that contain mirrors, and we conclude in Sec. |6} 

2 Absolute optical instruments 

In this section we recall some properties of absolute instruments from the point of view of 
geometrical optics, which will be useful for the subsequent calculations. We will consider 
radially symmetric AIs with the refractive index distribution n{r). In addition, we will focus 
on a specific class of absolute instruments, namely AIs of the first type [Ij. AI of the first type 
is a device with the property that every point A from its whole volume has a full image at 
some point B, which means that all rays emerging from A reach B. Since the role of the points 
A and B can be interchanged, it is clear that any ray emerging from the point A returns back 
to this point again, so A is an image of itself. In some AIs, it is the only image, in other ones 
there may be more images. 

The AI can be either three- or two-dimensional. In the latter case we consider a 2D prop- 
agation of rays in a 2D refractive index profile. Due to the radial symmetry of the device, the 
following quantity analogous to the mechanical angular momentum is conserved [1]: 

L = rn{r) sin a , (1) 

where a is the angle between the tangent to the particle trajectory and the radius vector. In 
3D, another consequence of the radial symmetry of the device is that each ray lies in a plane 
containing the centre of symmetry of the device O, so the ray trajectory can be described in 
polar coordinates (r, ip) both in the 2D and 3D cases. For nonzero L, the polar angle ip increases 
(or decreases) monotonically while r oscillates between the turning points r_ and r^, r^ < rj^. 
We will assume that for each possible value of L, there are just two turning points that coalesce 
when L reaches its largest possible value Lq. 

In order to ensure that any ray emerging from a point A eventually returns there, the change 
of the polar angle A(f corresponding to r changing from r_ to r+ must be a rational multiple 
of 27r, so we can write Aip = 7r//i, /i G (Q^ We will assume in this paper that the value of fx 
is the same for all possible angular momenta L, which is the most usual case. With respect to 
this assumption and since any point is an image of itself, the optical path ^aa = j ndl from A 
back to A is equal for all rays [T]. At the same time, the ray from A back to A consists of an 
integer number of segments on each of which r changes from r_ to r+ or back. Therefore also 
the optical path S between two turning points has to be equal for all rays, i.e., it is independent 
of the angular momentum L. This optical path can in general be expresses as 

SiL) = ri^= r ^dr 

J,_ cosa J,_ ^l - {L/nrf ^ ' 



^In [4j, m was used instead of /i. 



where we have used Eq. ([I| and written the dependence on L exphcitly. The fact that for 
absolute instruments S is independent of L is of a key importance for calculation of the spectra 
with the WKB method, as we will see in the following section. 

As the last thing we will express the optical path between the turning points in a different 
way. When L increases, the turning points r± approach each other and finally meet at the 
point To, the radius of the circular ray, for the maximum possible value of angular momentum 
L = Lq. The corresponding optical path S is then equal to the geometrical path A(fro = vrro//i 
multiplied by the refractive index n(ro), i.e.. 



S 



TTron{ro) 



(3) 



3 WKB calculation of spectra of AIs 

For simplicity we will consider a monochromatic scalar wave in the AI that can be described 
by the Helmholtz equation 

Aijj + k'^n{rf^lj = . (4) 

If the speed of light is set to unity, k is at the same time equal to the frequency u of the 
wave. Separating the radial and angular parts in terms of ilj{r,6,(f) = R{r)Yim{0,f) (in 3D) 
or ijj{r,ip) = i?(r)e^™''^ (in 2D) and making the substitution R{r) = r~^w{r) (in 3D) or R{r) = 
r~^/'^w{r) (in 2D), we get the equation for w{r) in the form 



w" + 



l{l + l) 



+ k^n 



2^2 



W 



(3D), 



w + 



m^ - 1/4 



+ k^n 



2„2 



w = (2D) , (5) 



where the prime denotes a derivative with respect to r. 

The equations ([s]) are analogous to the radial part of the Schrodinger equation obtained 
when solving a quantum-mechanical motion in a central potential. As was shown by Langer [H], 
however, a direct application of the WKB method to that equation leads to problems due to 
the centrifugal potential and termination of the r-axis at r = 0, and yields wrong results for the 
energy spectrum of the system. A careful analysis shows that the problem can be eliminated 
by the substitution /(/ + 1) — )■ (/ + 1/2)^ (or m? — 1/4 — )■ m? in 2D) — the so-called Langer 
modification, see [HI |15] or § 49 of |16j . Then the WKB method can be applied to the resulting 
equation directly and gives correct results. Exactly the same argument applies also in our case, 
which leads to the equations 



w" + 



{I + 1/2)^ 



+ k^n 



2^2 



W 



(3D) 



and 



w" + 



m 



TT + k^n 



2^2 



u; = (2D) 



(6) 



(7) 



We can now apply the WKB approximation to these equations directly. In the following we 
will write the formulas just for the 3D case as the treatment of the 2D case is analogous and 
can be obtained from the 3D results by the substitution / + 1/2 — > m. 
We first substitute w = pe^'^ in Eq. (p|, which yields 



+ i(2pV' + P^) 



,V_(l±ml 







(8) 



Then we separate the real and imaginary parts and in the real part neglect the term p" with 
respect to p0'^; this corresponds to the first order of WKB approximation. This way we get 
two equations 

0' = ±p(r), p0" + 2p'0' = O, (9) 

where we have denoted 

p{r) = JkM^ - ^^ ^ y^^' . (10) 

The first equation in rt9]) has the solution 0(r) = ±Jp{r)dr, the second one has the solution 
p = const X p~^'^. Combining these solutions, we get 

w{r) = — -1= exp I ±i / p{r) dr I . (11) 

vIpI \ J / 

The function w{r) has an oscillatory behaviour in the region where p is real, i.e., between the 
turning points r± at which p turns to zero. For r < r^ and r > r+, p is imaginary and the 
function w{r) describes evanescent waves in these regions. We now need to match the oscillatory 
and evanescent solutions in the three regions, which will give the quantisation condition for k. 
This matching, however, cannot be done within the WKB method itself because the condition 
\p"\ -C |p0'^| breaks up at the turning points. Still there are several ways how to get around 
this problem [16] . One of them is to treat r formally as a complex variable and circumvent the 
turning point in the complex r-plane. Another standard method is to approximate the bracket 
in Eq. d6|, i.e., p^, by a linear function of r in the vicinity of the turning point (say r+). Both 
of these methods give the same well-known result [H]: the wave decaying for r > r^ matches 
the following wave in the oscillatory region: 

const f r / i\ ^ I '^\ 

cos I / p(r ) dr H — . (12) 



This will be discussed in detail for a more general case in Sec. [5j A similar consideration can 
be made for the other turning point r_. To satisfy both the phase conditions simultaneously, 
it must hold 

^{l,k) = T:(N+]\ , Ar = 0,l,2,..., (13) 



2 
where we have defined 



$(/,fc)= / * p{r)dr= f * \k^n?- ^^^Y^^^ dr. (14) 






Eq. (13) is equivalent to the Born- Sommerf eld quantisation rule in quantum mechanics [Id] and 



together with Eq. ( 14 ) it presents the main result of this section. 

In some situations the condition imposed on the function w{r) may be different. In par- 
ticular, if there is a mirror in the AI at some radius r = 6, then the wave must vanish there. 



This condition will then lead to a different phase shift in Eq. (12) and hence to a shift in the 
eigenfrequencies, which will be discussed in Sec. |5j 



3.1 Calculation of the integral in Eq. (13) 



To find the eigenfrequencies for which the condition (13) is satisfied, we need to evaluate the 



integral (14). Quite remarkably, $(/,A;) is closely related to the optical path S in the (r, yj) 



plane between the turning points, defined by Eq. ^. It is easy to check that 

d<^{l,k) ^f 1 + 1/2' 



dk 



S 



(15) 



As we have seen in Sec. pi in AIs S{L) is independent of L, and therefore the right-hand side 
of Eq. (15) does not depend on k for a fixed /. This allows us to integrate Eq. (15) to get 



<!>{l,k) = Sk + c{l) 



(16) 



where c{l) is an integration constant with respect to k. To find it, we consider the situation 
when the turning points r± (given by the condition p = 0) coincide. As we have seen in Sec.j2| 
this corresponds to the value L = ron(ro) in Eq. ([2]). Comparing the square roots in Eqs. (|2| 

(/ + l/2)/[ron(ro)]. Moreover, in this 
-(/+l/2)^/(ron(ro)). 



and (10), we see that the corresponding value of k is ko 

because r^. = r_. Eq. (16) then yields the constant c(/) 



situation $ 

Substituting this to Eq. (16) and using Eq. (|3|, we get for k 



rori(ro) 



/i$ / 1 



(17) 



Combining Eqs. (17) and (13), we get finally the quantisation condition for the eigenfrequencies 



of a 3D absolute instrument 



k 



,{N+l)+l + l 



N 



rora(ro) 
The spectrum of a 2D AI is obtained by replacing / + 1/2 by \m\: 



0,1,2,... 



k 



ron(ro) 



fi N + -] +\m 



N = 0,l,2, 



;i8) 



(19) 



The formulas (18), (19) for the WKB-approximated spectrum of absolute instruments form the 



central result of this paper. They are based on the key relation ( 15 ) and on the fact that in 
AIs the optical path S is independent of the angular momentum. 

Inspection of the formulas (18), (19) shows two things. First, for absolute instruments 



the WKB spectrum is degenerate, and this degeneracy increases with the frequency. This is 
obvious from the fact that /x is a rational number for AIs, and usually it is a ratio of small 
integers. Then we can find a number of combinations (/, N) yielding the same frequency k, 



and this number increases with increasing / and/or A^. Since the formulas (18) and (19) give 
the WKB approximation to the spectrum, the exact spectrum exhibits this degeneracy only 
approximately and the eigenfrequencies form tight groups. The second property that follows 



from Eqs. (18) and (19) is that the frequency groups are positioned equidistantly because /x 
is rational and m,l and A^ all change in integer steps. As was shown in ^, both of these 
properties of the spectrum have key importance for the ability of AIs to focus waves. Eqs. ( 18 ) 



and ( 19 ) not only confirm the previous results, but in addition enable to calculate the absolute 



position of the spectral structure, which was not possible by the previous method. 



4 Examples 

To see how well the WKB spectrum approximates the actual spectrum, let us now apply our 
results on several examples of AIs. We will represent the spectra graphically by the function 



/c(z/) that has been introduced in [13] (the notation Ci;(z/) was used instead). Its value at 
integer v is simply the v^"^ eigenfrequency and for the non-integer z/, kiv) is given by linear 
interpolation |13]. Degenerate levels are clearly exhibited in the graph of kiv) as intervals 
where the function is constant. 

• Maxwell's fish eye, n{r 

WKB spectrum in 2D 



r) = 2/(1 + r^). In this case /i = 1, tq = 


= 1, which yields the 


A;wKB = iV + m + 1/2 . 


(20) 


spectrum is 




^analytic = \/ {^ + |m|)(A^ + |m + 1) . 


(21) 



With the exception of the lowest level, Eq. (20) approximates the correct values (21) very 



well: for A^ + \m\ = 1 the relative error is only 6%, for N + \m\ = 3 it is 1%, and it 
further decreases with increasing A^ + \m\. The two spectra are compared graphically in 
Fig. [I] (a). For the 3D Maxwell's fish eye profile, the WKB and exact eigenfrequencies 



are given by the same equations (20) and (21), respectively, with \m\ replaced by / + 1/2. 
Again the agreement is very good. 

Hooke index profile with the refractive index n(r) = \/2 — r^. This formula for the 
index is used also for r > \/2 where n becomes imaginary. This allows to calculate the 
eigenmodes and eigenfrequencies analytically as well; they are given by the Laguerre- 
Gaussian functions similarly as stationary states of an isotropic 2D harmonic oscillator 
in quantum mechanics. In this case yU = 2, tq = 1, which yields the 2D WKB spectrum 

kwKB = 2N+\m\ + l. (22) 

This is at the same time also the analytical form of the spectrum, so the WKB approxi- 
mation yields the exact values for the eigenfrequencies in this case. This is in a complete 
analogy with the situation in quantum mechanics where the WKB approximation also 
yields the exact values for the energy eigenvalues of a 2D harmonic oscillator. For the 3D 



Hooke profile both spectra are given by Eq. (22), again with \m\ replaced by / -|- 1/2; the 
agreement is perfect again. 



Kepler index profile, n{r) = y 2/r — 1. In this case /z = 1, tq = 1, which yields the 2D 



WKB spectrum in the form (20). This is at the same time also the analytical spectrum 
(again the formula for n{r) is used also for r > 2), so again the WKB gives here the exact 
spectrum. The spectra are shown in Fig. [I] (b). The situation in 3D is analogous. 

2D Minano lens, n(r) = 1 for r < 1 and n{r) = ^/2/r — 1 for r > 1. In this case 
/i = 2, To = 1, which yields the WKB spectrum in the form (22). Comparison with the 



numerically calculated spectrum is shown in Fig. [I] (c). The agreement is very good 

• As the last example we consider a 2D AI with the refractive index n{r) = 2r^^{[(2 — 
r)/r]'^/^ _l_ [7-/(2 — r)]'^/^}~^ derived from the general formula (18) of [4J with /(r) = 2 — r 
and /i = 3. This yields for the WKB spectrum 

kwKB = SN+\m\ + -, (23) 

which again approximates the exact spectrum well, as can be seen in Fig. [I] (d). 

We see that in all these cases the agreement of the WKB spectrum with the exact one is very 
good, even for the lowest states. This is quite a remarkable feature of the WKB approximation 
that our optical situation shares with quantum mechanics. 
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Figure 1: Comparison of the spectra of several AIs calculated with the WKB method (blue 
line) with the analytically or numerically calculated spectra (red dashed line). The spectra are 
represented by the function fc(z/), degenerate levels correspond to plateaus in the graph, (a) 
Maxwell's fish eye, (b) Kepler profile, (c) Minano lens, (d) a 2D AI with the refractive index 
n(r) = 2r~^{[(2 — r)/r]^/'^ + [r/(2 — r)]^/^}"^. In (b), the analytical and WKB spectra are 
identical. 



5 Mirrors 

So far we have considered absolute instruments without mirrors, so we required the wave in the 
evanescent region to decay gradually. However, in many AIs there are mirrors that constitute 
their important parts. The mirrors play one or both of the following two roles. First, in some 
AIs they limit the size of the device, e.g. from infinite size to a finite one. This is the case of 
Maxwell's fish eye mirror (MFEM) [Hill]. The second possible role is that mirrors can eliminate 
regions where the refractive index goes to zero; this is again the case of MFEM or e.g. of the 
modified Minano lens |1]. 

Consider a situation that there is a spherical (or in 2D, circular) mirror at radius r = 6 in 
the AI. Due to the boundary condition the wave has to vanish at it, i.e., w{b) = 0. We have to 
distinguish two cases: either the mirror is placed in the evanescent or in the oscillatory region. 
In the following we will treat them separately. 

5.1 Mirror in the evanescent region 

Suppose first that the mirror is placed in the evanescent region, e.g. fe > r+. Then for r > r+ 



we need not only the exponentially decaying solution in Eq. (11) but also the exponentially 
growing one to create their superposition that vanishes at r = 6. This will cause a phase shift 
of the wave in the oscillatory region and consequently a shift in the spectrum compared to the 
situation without the mirror. 

To describe such a situation, we will use one of the standard methods for matching the WKB 
solutions in the oscillatory and evanescent regions employing Airy functions. For this purpose 
we approximate the bracket in Eq. (^, i.e., p"^, by a linear function of r in the vicinity of the 
turning points. To write this explicitly for the turning point r_|_, we put p^(r) ^ —a{r — r+), 
for which Eq. (ml) has exact solutions Wi{r) = Ai(a^'^(r — r+)) and W2{r) = Bi(a^'^(r — r+)) 
(up to a multiplicative constant). Using the asymptotic formulas for Airy functions [IT] that 
are very accurate for |x| > 2, namely 

^ , -x)-V4 sin[i(-x)3/^ + vr/4] (x < 0) 




{x > 0) 



and 



n 



-l/2r_^^-l/4p^^r2/_^A3/2 



x)-i/^cos[f(-x)3/^ + 7r/4] (x<0) 



Bi{x)^i _,'_,, ,2 .L ^ . (25) 



n 



-i/2a;-i/4exp(|x3/2) (x > 0) 



we can express the solutions 14^1,2 (^) for r < r+ as 

,1/6 



Wi{r) ^ cos I / p{r')dr'-\ — ) (26) 

V^ \Jr+ 4 J 

fpir')dr' + ^^ , (27) 



aV6 ^ / r 



W2[r) ~ sm 

y/Tip 



where the linearisation of p^ was used again to write the expressions (ix)'^/^ in Eqs. (24) 



and (25) as integrals of p. Obviously, the solutions (26) and (27) can be expressed as linear 
combinations of the solution s ([TT] ) in the oscillatory region. This way Eqs. (26) and (27) describe 
the same solutions as Eq. (11), and a similar consideration can be made for the evanescent 



region. This allows us to express the approximate solutions of Eq. mh in the evanescent and 





w (r) 




(b) 



Figure 2: Example of the functions relevant for describing the approximate solutions of Eq. (|6|. 
(a) The function p^(r) chosen to be significantly non- linear, (b) the corresponding functions 
/i(r) (red) and /2(^) (dashed blue) from Eqs. (28). In (c) we compare a numerical solution of 
Eq. do]) satisfying w{2) = 0, w'{2) = —0.02 (red curve) with its approximate solution in terms of 
a superposition of the functions /i(r), f2{r) (dashed blue curve) satisfying the same conditions. 
The agreement between the two curves is very good, in particular the position of the nodes 
that has the key importance for the spectrum. The y-axis of the graphs is placed at the turning 
point r+ (in this case r+ = 1.3). 




f.ir) 





Figure 3: As in Fig. pi but with an even stronger non-linearity of the function p(r)^, and the 
parameters r^ = 1.5, w(0.5) = l,w'{0.5) = —1. 



oscillatory regions as well as in the vicinity of the turning point r+ by a single formula using the 



Airy functions Ai and Bi. Their argument is obtained by expressing the x variable in Eqs. (24) 



and (25) in terms of the integral of p with the help of Eqs. (26) and (27). This way we arrive 



at the following approximate solutions of Eq. ^ : 



/i(r)=Ai 




/2(r)=Bi 




(28) 



In this expression the integral is real or purely imaginary for r < r_(_ or r > r_|_, respectively. 
Its square in these regions is positive or negative, respectively, and so is the third root. It is 



easy to check that if p (r) = —a{r — r+), then the brackets in Eq. (28) turn into a 



.1/3/ 



for all r, so in that case Eqs. (28) represent the exact solutions. Moreover, even if p^(r) is 



non-linear, Eqs. (28) still provide a very good approximation to the solutions of Eq. ([6]). Most 
importantly, the positions of the nodes match very well the positions of the nodes of the exact 
solution. This can be seen on the examples in Figs. [2] and |3} 

Eqs. (28 ) enable to find the phase of the wave in the oscillatory region very accurately. What 



we have to do is to find the coefficients ci,2 such that the superposition w{r) = C1/1+C2/2 satisfies 
the condition w{b) = 0. The solution is, up to a common multiphcative factor, ci = Bi(,^) and 



C2 = — Ai(^), where ^ = ~ y {2 K^ P ^''") ■ Using the asymptotic forms (26), (27) in the 
oscillatory region, we find that in this region 



w{r) ^ — -^ cos I / p{r') dr' + rj] , (29) 



const 

Vp 



with the phase shift 



Bi(n - Ai(n , , 

^ = "^^'""Bi(0+Ai(0 ^'°^ 



The difference between r] and n/A gives a correction factor for Eqs. (18) and (19) that has to 
be added to the curly parentheses. A similar consideration can be made in the situation when 
there is a mirror in the other evanescent region r < r_. Naturally, the closer the turning point 
is to the mirror, the larger will be the correction. This way the largest shifts will be exhibited 
by the levels corresponding to the smallest values of m. 

Fig. |4] shows the results of our method on two examples of AIs where the mirror has been 
placed at the boundary of the classically accessible region at r = 6 = 2. The first one is 
the Kepler profile, the second one is Minano lens. The WKB spectrum is compared with the 
numerically calculated spectrum. The difference is so small that it can be hardly seen in the 
pictures. For comparison, the WKB spectrum in the absence of the mirror is also shown. The 
mirror obviously shifts all the levels up as could be expected. A closer inspection of the spectrum 
reveals that the shift indeed decreases with increasing m due to the increasing separation of 
the turning point from the mirror. 

Fig. |4] also clearly reveals that the spectrum is no more so degenerate as it was in the 
absence of the mirror. As was shown in |1], the regularity and degeneracy of the spectrum is 
connected with the quality of imaging by the device, so adding the mirror degrades somewhat 
the imaging, in particular at low frequencies. Note that this happens even though the mirror 
is situated beyond the turning point r+ for most rays (i.e., for all rays with the exception of 
the ones with L = for which the mirror is just at the turning point). In terms of geometrical 
optics, nothing changes when the mirror is added because only the rays with L = reach it, 
and even their trajectory remains unchanged. However, unlike the rays, the wave "feels" the 
mirror even if it is located in the evanescent, i.e., classically inaccessible region. 

5.2 Mirror in the oscillatory region 

We will now consider the situation when the mirror is not in the evanescent, but rather in the 
oscillatory region, which is the case e.g. of Maxwell's fish eye mirror. Suppose that the optical 
medium is inside the mirror, again as in the case of MFEM. The upper turning point r+ is 
then located at the mirror itself, so w{r) must vanish there and so for r < r^ the wave must 



behave according to Eq. (29) with rj = tx 12. This means an additional constant phase factor of 



7r/4 compared to Eq. (12), which then has to be added into the curly parentheses in Eqs. (18) 



and (19). Applying this resuh to the 2D MFEM (for which n{r) = 2/{l + r^), /i = 2 and b = 1] 



we get for the WKB spectrum 

fcwKB = 2N+ \m\ + 3/2 . (31) 

The analytical form of the spectrum is 

^analytic = V(2A^ + l"^l + 1) (2A^ + 1^1 + 2) ■ (32) 
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Figure 4: Comparison of the WKB spectra (blue) and numerically calculated spectra (dashed 
red) for two AIs with mirrors, (a) Kepler index profile surrounded by a mirror at radius 6 = 2. 
(b) Minano profile surrounded by a mirror at radius 6 = 2. The difference between the curves 
can be hardly seen, the relative errors make a few per cent even for the lowest levels and becomes 
completely negligible for higher frequencies. The thin black lines show the WKB spectra of the 
corresponding devices without the mirror. 

This way, the WKB method gives here very good results in a similar way as for the simple 
Maxwell's fish eye. 

6 Conclusions 

In this paper we have developed a method for calculating frequency spectra of absolute in- 



struments using the WKB approximation. Our method is based on a the key relation (15) 



between the optical path S and the semiclassical phase change $ between the turning points, 
and on the fact that in AIs of the first type the optical path S is independent of the angular 
momentum. The method proved to be very efficient for describing the spectrum accurately even 
for the lowest levels. The results confirmed the previously derived properties of the spectra, 
in particular that the frequencies are strongly degenerate in AIs and that they form almost 
equidistantly spaced groups. Since these properties of the spectrum have key importance for 
focusing of waves by AIs [13j, the WKB method shows that sharp focusing of rays in fact 
implies that waves will be focused as well. This way the WKB method provides a nice and 
important bridge between geometrical and wave optics of absolute instruments. 

We applied our method also to AIs that contain mirrors; there, the Airy functions turned 
out to be very helpful in describing the wave in both the oscillatory and evanescent regions, 
and to calculate the eigenfrequency shifts caused by the mirror. When the mirror is added into 
the evanescent region, it turns out that even though light rays are not infiuenced at all, the 
spectrum can be infiuenced strongly. In particular, the spectrum becomes less regular, which 
can somewhat degrade imaging by the absolute instrument. On the other hand, placing the 
mirror in the oscillatory region leads to a constant shift in the spectrum. 

It may be somewhat surprising, but certainly very satisfactory, how accurate results the 
WKB method gives for the spectra of absolute instruments. The situation is similar to quantum 
mechanics where the WKB method also gives very precise values for the energy spectrum in 



11 



many situations. An interesting question for further research could be how the WKB method 
can be apphed to AIs of the second type where not all points have their full images, or to even 
more general devices. 
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